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Abstract. We explore the cosmological implications of five modified gravity (MG) models 
by using the recent cosmological observational data, including the recently released SNLS3 
type la supernovae sample, the cosmic microwave background anisotropy data from the 
Wilkinson Microwave Anisotropy Probe 7-yr observations, the baryon acoustic oscillation 
results from the Sloan Digital Sky Survey data release 7, and the latest Hubble constant 
measurement utilizing the Wide Field Camera 3 on the Hubble Space Telescope. The MG 
models considered include the Dvali-Gabadadze-Porrati(DGP) model, two f(R) models, and 
two f(T) models. We find that compared with the ACDM model, MG models can not lead 
to a appreciable reduction of the Xmin- The analysis of AIC and BIC shows that the simplest 
cosmological constant model(ACDM) is still most preferred by the current data, and the 
DGP model is strongly disfavored. In addition, from the observational constraints, we also 
reconstruct the evolutions of the growth factor in these models. We find that the current 
available growth factor data are not enough to distinguish these MG models from the ACDM 
model. 
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1 Introduction 

Since its discovery in 1998 [1], the cosmic acceleration has become one of the central prob- 
lems in theoretical physics and cosmology. The most common approach to explain the cosmic 
acceleration is modifying the right hand side of Einstein equation, i.e., introducing a myste- 
rious composition called dark energy (DE) [2]. Numerous theoretical DE models have been 
proposed in the last decade [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13], while it is still as daunting as 
ever to understand the nature of the cosmic acceleration [14]. 

In addition to the introduction of the DE component, another popular route to explain 
the cosmic acceleration, i.e., modifying the left hand side of Einstein equation, has also drawn 
a lot of attentions. They are the so-called modified gravity (MG) models [15, 16]. So far, a lot 
of MG models have been proposed, providing various interesting mechanisms for the cosmic 
acceleration. For example, in the DGP braneworld model [17], gravity is altered at immense 
distances by slow leakage of gravity off from our three-dimensional universe, leading to the 
apparent acceleration in the cosmological scale. In the f{R) gravity [18-21], the Ricci scalar 
R is generalized to a function f(R), providing another possible mechanism for the cosmic 
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acceleration. A similar idea is the f(T) gravity, where the torsion T is replaced by a function 
f(T). Other popular MG models include MOND [22], TeVeS [23], Brans-Dicke gravity [24], 
Gauss-Bonnet gravity [25], Horava-Lifshitz gravity [26], and so on. For recent reviews about 
MG models, see e.g. [16, 27]. 

There are mainly two ways to constrain the cosmological models by using the observa- 
tional data. First, they can be tested through their effects in the cosmic expansion history, 
which can be measured through Type la supernoave (SNIa), baryon acoustic oscillations 
(BAO), and so on. Second, they can be tested through their effects in the cosmic growth 
history, which can be measured through the Integrated Sachs Wolfe (ISW) effect [28], weak 
lensing, galaxy clusters, and so on. So far, many studies have been performed [29-52]. 

In this work, we will mainly test the MG models against the observations of the cosmic 
expansion history. Recently, a high-quality joint sample of 472 supernovae (SN), the SNLS3 
SNIa dataset [53], was released. Moreover, the systematic uncertainties of the SNIa data were 
nicely handled in work [53]. It is interesting to study the performance of MG models when 
confronted with this newly released SNIa dataset. Thus, in this paper we will adopt SNLS3 
dataset and systematically study 5 representative MG models, including DGP model, two 
f(R) models, and two f(T) models. To perform a comprehensive analysis, we also include 
the cosmic microwave background (CMB) anisotropy data from the Wilkinson Microwave 
Anisotropy Probe 7-yr (WMAP7) observations [59], the baryon acoustic oscillation (BAO) 
results from the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) [60], and the latest 
Hubble constant measurement from the Wide Field Camera 3 (WFC3) on the HST [61]. 
Utilizing these data, we constrain the parameter space of these MG models, and compare 
them with the ACDM model by using the AIC [62] and BIC [63] criterion. 

An important tool to study the cosmic growth history is growth factor [64]. However, 
compared with the cosmic expansion history data, the current growth factor observations 
have larger errors (their redshifts are also not well determined). Moreover, some data points 
of growth factor observations are obtained by assuming ACDM, and therefore their use to 
test models different from ACDM may not be reliable. So in this work we will not include the 
growth factor data when performing observational constraints on the MG models. Instead, 
making use of the constraints on the MG models obtained from the observations of the cosmic 
expansion history, we reconstruct the evolutions of the growth factor along with the redshift 
z in these models and compare the results with the growth factor data. In this way, we can 
see whether it is possible to distinguish these MG models from the growth factor data. 

This paper is organized as follows: In Sec. II, we give a brief and overall description on 
the models considered here, as well as the analysis methods used. In Sec. Ill, we introduce in 
detail the current observational data of the cosmic expansion history and the cosmic growth 
history, respectively. In Sec. IV, we study the cosmological interpretations of the MG models 
considered here and give a summary of these models by making a comparison among them. 
In Sec. V, we give a short discussion. In addition, we present two appendixes about the 
initial condition of evolution equation of growth factor and the effects of the parameter n in 
f(R) model. In this work, we assume today's scale factor oo = 1, so the redshift z = a" 1 — 1; 
the subscript "0" always indicates the present value of the corresponding quantity, and the 
unit with c = H = 1 is used. 
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2 Models and Methodology 



2.1 Models 

In this work, we will investigate 5 representative MG models, including 

1. The DGP model 

2. A power-law type f(T) model {f(T)pi) 

3. An exponential type f(T) model (/(T)exp) 

4. The f(R) model proposed by Starobinsky (f(R)st) 

5. The f(R) model proposed by Hu and Sawicki (f(R)iis) 

The DGP model is a braneworld model [17], where gravity is altered at immense dis- 
tances by slow leakage of gravity off from our three-dimensional universe. In the f(T) models, 
the torsion T in the Lagrangian density is replaced by a generalized function T + f(T). In 
this work, we consider two f(T) models, in which the functions f(T) are of power-law type 
[20] (hereafter /(T)pl model) and exponential type [21] (hereafter /(T)exp model), respec- 
tively. Similarly, in the f(R) models, the Ricci scalar R is replaced by R + f(R), and we 
will investigate the two popular f(R) models proposed by Starobinsky [19] (hereafter f(R)st 
model) and by Hu and Sawicki [18] (hereafter f(R)HS model), respectively. We will give a 
detailed introduction about the formulas of these models in Sec. IV. 

2.2 The x 2 analysis 

To study the cosmological interpretations of the models listed in the previous subsection, 
we have to make use of the data from the recent cosmological observations. In practice, 
theoretical models and observational data can be related through the x 2 statistic. For a 
physical quantity £ with experimentally measured value £ b s , standard deviation a^ and 
theoretically predicted value ^th, the x 2 takes the form 

Xf = -2 • (2- 1 ) 

In case that there are more than one independent physical quantities, one can construct total 
X 2 by simply summing up all the x| s , i- e -> 

X 2 = E4" ( 2 - 2 ) 

t 

In the following subsections, we will present the detailed forms of the x 2 function of the 
cosmological observations considered in this work. 

Once having got the x 2 function, one can determine not only the best-fit model pa- 
rameters but also the la and the 2a confidence level (CL) ranges of the model considered 
by minimizing the total x 2 an d calculating Ax 2 = X 2 ~ X 2 m ini respectively. This need a 
thorough exploration of the parameter space of the x 2 function. This procedure is commonly 
accomplished through the Monte Carlo Markov chain (MCMC) technique. In this work, we 
make use of the publicly available CosmoMC package [65] and generate O(10 6 ) samples for 
each set of results presented in this paper. 

2.3 Model Comparison 

A statistical variable must be chosen to enforce a comparison between different models. The 
Xmin may be the simplest one, but it is difficult to compare models with different number of 
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parameters. Thus, in this work we use the commonly used AIC and BIC as model selection 
criterions. 

The information criteria (IC) is the most frequently used method to assess different 
models. It is also based on a likelihood method. In this paper, we use the Akaike information 
criterion (AIC) [62] and the Bayesian information criterion (BIC) [63] as model selection 
criteria. These criteria favor models that have fewer parameters while giving a same fit, and 
have been applied frequently in the cosmological contexts [66-69]. 

The AIC [62] is defined as 

AIC = —2hiC max + 2n p , (2.3) 

here C max is the maximum likelihood and n p denote the number of free model parameters. 
Note that for Gaussian errors, C max is related to the Xmin by 

-2lnC max = X 2 min , (2-4) 

so the difference in AIC can be simplified to AAIC = AxLin + 2An p . As mentioned in Ref. 
[71], there is a version of the AIC corrected for small sample sizes, AIC C = AIC + 2n p (n p — 
1)/(N — n p — 1), which is important for N/n p < 40. In our case, this correction is negligible. 
The BIC, also known as the Schwarz information criterion [63], is given by 

BIC = -2 In C max + n p In N, (2.5) 

where N is the number of data. Also, in this case, the absolute value of the criterion is not 
of interest, only the relative value between different models, ABIC = Ax^ n + An p lnN, is 
useful, it should be mentioned that the AIC is more lenient than BIC on models with extra 
parameters for any likely data set IniV > 2. Generally speaking, a difference of 2 in BIC 
(ABIC) is considered as positive evidence against the model with the higher BIC, while a 
ABIC of 6 is considered as strong evidence. 

It should be noted that the IC alone can only imply that a more complex model is not 
necessary to explain the data, since a poor information criterion result might arise from the 
fact that the current data are too limited to constrain the extra parameters in this complex 
model, and it might become preferred with the future data. Thus, this only reflects the 
current situation for the theoretical models. 

3 Cosmological Observations 

In this section, firstly, we will introduce in detail the observational data of the cosmic expan- 
sion history, and then we will present the observational data of the cosmic growth history. 

3.1 Observations of the Cosmic Expansion History 

In this paper, we will constrain the MG models by using the observational data of the 
cosmic expansion history, including the recently released SNLS3 SNIa sample [53] , the CMB 
anisotropy data from the WMAP7 observations [59], the BAO results from the SDSS DR7 
[60], and the latest Hubble constant measurement utilizing the WFC3 on the HST [61]. In 
the following, we will include these observations into the x 2 analysis. 



-4- 



3.1.1 SNIa 



A most powerful probe of the cosmic acceleration is Type la supernovae (SNIa), which are 
used as cosmological standard candles to directly measure the cosmic expansion. Recently, a 
high-quality joint sample of 472 supernovae (SN), the SNLS3 SNIa dataset [53], was released. 
This sample includes 242 SN at 0.08 < z < 1.06 from the Supernova Legacy Survey (SNLS) 3- 
yr observations [54], 123 SN at low redshifts [55] [56], 93 SN at intermediate redshifts from the 
Sloan Digital Sky Survey (SDSS)-II SN search [57], and 14 SN at z > 0.8 from Hubble Space 
Telescope (HST) [58]. This SNIa sample has been used to probe properties of DE [72] [73] [74] , 
to constrain the cosmological parameters [75], and to test the cosmological models [76]. Here, 
we briefly discuss the x 2 function of the SNLS3 SNIa data, which can be downloaded from 
[77]. 

The x 2 function of the SNLS3 SNIa dataset takes the form 

x| N = A^ T • Cr 1 • (3.1) 

where C is a 472 x 472 covariance matrix capturing the statistic and systematic uncertainties 
of the SNIa sample, and Arit = rn^ — m mo d is a vector of model residuals of the SNIa sample. 
Here is the rest-frame peak B band magnitude of the SNIa, and m mo d is the predicted 
magnitude of the SNIa given by the cosmological model and two other quantities (stretch 
and color) describing the light-curve of the particular SNIa. The model magnitude m mo d is 
given by 

mmod = 5log w V L (z hc i,z cmb ) - a(s - 1) + PC + M. (3.2) 

Here T>l is the Hubble-constant free luminosity distance. In a spatially flat Friedmann- 
Robertson- Walker (FRW) universe (the assumption of flatness is motivated by the inflation 
scenario), it takes the form 

T>L(z heh Z zcmh ) = (1 + Zhel) / J-, / ,, ■ (3.3) 

Jo ) 

Here £ cm b and Zhel are the CMB frame and heliocentric redshifts of the SN, s is the stretch 
measure for the SN, and C is the color measure for the SN. E{z') = H(z)/Hq, and H{z) is 
the Hubble parameter. Notice that different model will give different H(z) and E{z). a and 
/3 are nuisance parameters which characterize the stretch-luminosity and color-luminosity 
relationships, respectively. Following [53], we treat a and f3 as free parameters and let them 
run freely. 

The quantity M in Eq. (3.2) is a nuisance parameter representing some combination 
of the absolute magnitude of a fiducial SNIa and the Hubble constant. In this work, we 
marginalize A4 following the complicated formula in Appendix C of [53]. This procedure 
includes the host-galaxy information [78] in the cosmological fits by splitting the samples 
into two parts and allowing the absolute magnitude to be different between these two parts. 

The total covariance matrix C in Eq. (3.2) captures both the statistical and systematic 
uncertainties of the SNIa data. One can decompose it as [53], 

C = D s tat + C s tat + C sys , (3.4) 

where D s tat is the purely diagonal part of the statistical uncertainties, C s t a t is the off-diagonal 
part of the statistical uncertainties, and C sys is the part capturing the systematic uncertain- 
ties. It should be mentioned that, for different a and (3, these covariance matrices are also 
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different. Therefore, in practice one has to reconstruct the covariance matrix C for the cor- 
responding values of a and /3, and calculate its inversion. For simplicity, we do not describe 
these covariance matrices one by one. One can refer to the original paper [53] and the public 
code [77] for more details about the explicit forms of the covariance matrices and the details 
of the calculation of x| N . 

3.1.2 CMB 

Then we turn to the CMB observational data. We use the " WMAP distance priors" given by 
the 7-yr WMAP observations [59]. The distance priors include the "acoustic scale" I a, the 
"shift parameter" R, and the redshift of the decoupling epoch of photons . The acoustic 
scale I a, which represents the CMB multipole corresponding to the location of the acoustic 
peak, is defined as [59] 

Ia={1 + z*) — r— . (3.5) 

r s {z*) 

Here Da{z) is the proper angular diameter distance, given by 

„ , x 1 f z dz' 

and r s (z) is the comoving sound horizon size, given by 

, , 1 f 1/(1+z) da 

r s (zj = —p= I :, (3.7) 

VsJo a 2 if(o)Vl + (3n6o/4n 70 )o 

where fi&o and fi 7 o are the present baryon and photon density parameters, respectively. In 
this paper, we adopt the best-fit values, Q 60 = 0.02253/i" 2 and ft 7 o = 2.469 x W~ 5 h~ 2 (for 
T C mb = 2.725 K), given by the 7-yr WMAP observations [59]. Here h is the reduced Hubble 
parameter satisfying Hq = h ■ 100 km/s/Mpc. The fitting function of z* was proposed by Hu 
and Sugiyama [79]: 

z* = 1048[1 +0.00124(^ /i 2 )-°- 738 ][l + gi(n m0 h 2 ) 92 ], (3.8) 

where 

_ 0.0783(^ /i 2 )- - 238 0-560 
91 l + 39.5{n b0 h 2 ) ■ 7 ™ , 92 l + 21.1(fi 60 /j 2 ) 1 - 81 ' 1 1 

In addition, the shift parameter R is defined as [80] 



R(z*) = ^n m0 H${l + z*)D A (z*), (3.10) 

here £l m o is the present fractional matter density. 
As shown in [59], the x 2 °f the CMB data is 

Xcmb = ( x i bs ~ x i' l )(Cc , MB)«i( x j feS _ ^i' 1 )' (3-11) 

where Xi = (Ia,R,z*) is a vector, and (C^ MB )ij is the inverse covariance matrix. The 7- 
yr WMAP observations [59] had given the maximum likelihood values: Ia(z*) = 302.09, 
R(z*) = 1.725, and z* = 1091.3. The inverse covariance matrix was also given in [59] 

2.305 29.698 -1.333 
(C~] /B ) = | 29.698 6825.27 -113.180 | . (3.12) 
-1.333 -113.180 3.414 
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3.1.3 BAO 

Next we consider the BAO observational data. We use the distance measures from the SDSS 
DR7 [60]. One effective distance measure is the Dy(z), which can be obtained from the 
spherical average [81] 

1/3 

(3.13) 



D v {z) = 



{l + zYDi(z)- 



'H(z) 

where Da (z) is the proper angular diameter distance. In this work we use two quantities 
do.2 = r s {zd)/ Dy (0.2) and c?o.35 = r s {z ( j) / 'Dy(0.35). The expression of r s is given in Eq.(3.7), 
and Zd denotes the redshift of the drag epoch, whose fitting formula is proposed by Eisenstein 
and Hu [82] 

r 

i + h(n b0 h 2 ) b2 , (3.14) 



1291(fi m0 /i 2 ) - 251 



1 + 0.659(O m0 /i 2 ) 



0.828 



where 



h = 0.313(a m0 /i 2 )" - 419 [1 + 0.607(O m0 /i 2 ) ' 674 ] , (3.15) 

b 2 = 0.238(n m0 h 2 ) - 223 . (3.16) 

Following [60], we write the x 2 for the BAO data as, 

Xbao = A Pi( C BAo)ii A Pji ( 3 - 17 ) 

where 

A Pi = pf ata - Pi , pf ta = < a 2 ta = 0.1905, pf t& = < a 3* 5 a = 0.1097, (3.18) 
and the inverse covariance matrix takes the form 

! / 30124 -17227\ 

\ BAO> \ T70ov sfin77 • {o.iy) 



'BAO> ^-17227 86977 J 
3.1.4 The Hubble Constant H 

We also use the prior on the Hubble constant Hq. The precise measurements of Hq will be 
helpful to break the degeneracy between it and the DE parameters [83]. When combined 
with the CMB measurement, it can lead to precise measurement of the DE EOS w [84]. 
Recently, using the WFC3 on the HST, Riess et al. obtained an accurate determination of 
the Hubble constant [61] 

H = 73.8 ± 2.4km/s/Mpc, (3.20) 
corresponding to a 3.3% uncertainty. So the x 2 °f the Hubble constant data is 



h - 0.738 



2 



^=( :l ao2f : J' < 321) 

3.1.5 Combining of the SNIa, CMB, BAO and Hq data 

Since the SNIa, CMB, BAO and Hq are effectively independent measurements, we can com- 
bine them by simply adding together the x 2 functions, i.e., 

Xau = Xsn + Xc'mb + Xbao + xl- (3-22) 
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3.2 Observations of the Cosmic Growth History 

As mentioned above, growth factor is an important tool to study the cosmic growth history. 
In the following, we will give a brief introduction to the growth factor. 

As shown in [85], in the subhorizon limit, the evolution equation of the matter density 
perturbation of MG theory has the following form, 



5 + 2H5 - ATiG effPm 5 = 0, (3.23) 



where 5 = ^ SL , p m is the matter energy density and G e ff is the effective Newton's constant. 
The explicit forms of G e ff for the MG models considered in this paper will be given in Sec. 4. 
The growth factor / is defined as 

/ - (3.24) 
a ma 

Thus, Eq.(3.23) becomes 

* + / 2 + (^ + 2)/ = |5^. (3.25, 



din a K H ,J 2 3H 2 

To solve Eq.(3.25) numerically, the initial condition is taken as f(z = 30) = 1 (see Appendix 
A). Notice that different models correspond to different H{z) and G e ff, and thus give different 
f(z). We will make use of this equation to determine the predicted evolution history of the 
growth factor in MG models. 

To test the theoretical prediction of growth factor obtained from Eq.(3.25), it may be 
helpful to consider the current growth factor data [86-88], which are summarized and listed 
in the Table 1. Compared with the cosmic expansion history data, the current growth factor 
observations have larger errors, and their redshifts are also not well determined. Moreover, 
as pointed out in [87], some data points of Table 1 are obtained by assuming ACDM when 
converting redshifts to distances for the power spectra, and therefore it may be unreliable 
to use them in the models different from ACDM. So in this work, we will not include the 
growth factor data in our y 2 analysis. Instead, by using the constraints obtained from the 
cosmic expansion history data, we will explore the evolutions of these models' growth factor 
along with redshift z, and compare the results with the growth factor data. 

4 Results and Discussions 

In the following, we firstly constrain the parameter spaces of the models by using the observa- 
tional data of the cosmic expansion history, then we plot the corresponding evolutions of these 
models' growth factor along with redshift z, and compare the results with the observational 
data of growth factor f f, s . 

4.1 DGP model 

We start with the flat DGP model [17]. In this model, the cosmic expansion history E(z) = 
H(z)/Hq is determined by the following equation [96]: 



e(z) = [^n m0 (i + zf + o r0 (i + zY + n rc + vo^J . (4.i) 

Here, Sl rc = (1 — Q, m o — ^ro) 2 /4, and Q r o is the present fractional radiation density given by 
[59] 

U r0 = 7 o(l + 0.2271iV e// ), ft 7 o = 2.469 x l(T 5 /i~ 2 , N eff = 3.04, (4.2) 
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z 


fobs 


References 


0.15 


0.49 ±0.1 


[88, 891 

L " 1 


0.35 


0.7 ±0.18 


[901 

L 1 


0.55 


0.75 ±0.18 


[911 

L J 


0.77 


0.91 ±0.36 


[881 

L J 


1.4 


0.9 ±0.24 


[921 

L J 


3.0 


1.46 ±0.29 


[93] 


2.125 - 2.72 


0.74 ±0.24 


[94] 


2.2 - 3 


0.99 ± 1.16 


[95] 


2.4-3.2 


1.13 ± 1.07 


[95] 


2.6 - 3.4 


1.66 ± 1.35 


[95] 


2.8 - 3.6 


1.43 ±1.34 


[95] 


3-3.8 


1.3 ±1.5 


[95] 



Table 1. The combined observational data on the growth factor /. 




Figure 1. LeftPanel: Probability contours at the 68.3% and 95.4% confidence levels in £l m0 ~h plane, 
for the DGP model. This figure shows that there is an inconsistency between different cosmological 
probes for DGP model. Right Panel: The evolutions of f(z) predicted by the ACDM and DGP 
models. The best-fit values and 2a errors regions determined by SNIa+CMB+BAO+i?o analysis are 
shown. 



where s2 7 o is the present fractional photon density, and N e jj is the effective number of 
neutrino species. 

In this model, the effective Newton constant G e ff is no longer the simple constant G, 
instead, it takes the following form [97], 

2(1 + 2^)) 

G *" = ° 3(1 + 0^))- (43) 

Like the ACDM model, DGP is also an one-parameter model with the sole parameter 
O m o- However, it has been shown to be disfavored by the cosmological observations [29, 30, 
98]. In this work we obtain similar result. In the left panel of Fig.l, we plot the contours 
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of 68.3% and 95.4% confidence levels in the £l m Q-h plane, for the DGP model. Constraints 
from SNLS3, CMB+H , BAO+if , and SNLS3+CMB+BAO+# are shown in contours 
with different colors. The figure shows that there is an inconsistency between different 
cosmological probes in the DGP model. This also implies that the DGP model is disfavored 
by these cosmological probes. At the 68.3% and 95.4% confidence levels, we get 

n — n 97c;q+0.0151 +0.0312 . _ n fi4A r+0.0093 +0.0189 (a a\ 

S'mO — U - Z '°"J-0.0140 -0.0274' n ~ u -°^°-0.0093 -0.0186- 

In the right panel of Fig.l, we plot the evolutions of f(z) predicted by the ACDM model 
and the DGP model along with the observational data of growth factor f j, s . From this figure, 
we can see clear difference between the predicted evolutions of f{z) from the DGP model and 
the ACDM model. However, the current available growth factor data are still not powerful 
enough to distinguish these two models. 

4.2 f(T) models 

Then, we turn to the f(T) models. The action of f(T) models is 

S = ^ j d^x^-g [T + f(T)] + S m + S r , (4.5) 

here, S m and S r are the action of the matter content and the radiation content, and T is the 
torsion scalar [99]. 

Assuming a flat homogeneous and isotropic FRW universe, T satisfies 

T = 6H 2 , (4.6) 

and the modified Friedmann equations are obtained as [21] 

H* = ^( Pm + Pr )-l-2Hif T , (4.7) 

3 6 

• l W 2 + f + l2H*f TT 

4 1 + f T - I2m f TT ' 

here and throughout, p m and p r are matter density and radiation density, fx and /tt are 
defined as 

* = hr^%. (4.9) 
For f(T) models, G e ff is given by [100] 

G *ff = 7TT- 

1 J JT 

Notice that, if f(T) is a constant, the term f(T) acts just like a cosmological constant. 

In this paper, we consider two f(T) models with different types of parametrization: One 
is the power law model considered in [20], the other is exponential form model proposed by 
Linder [21]. For simplicity, hereafter we will call them f(T)pi, and f(T)EXP, respectively. 
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Figure 2. Left Panel: Probability contours at the 68.3% and 95.4% confidence levels in fi m o-n 
plane, for the /(T)pl model. Right Panel: The evolutions of f(z) predicted by the ACDM and 
f{T)pL models. The best-fit values and 2a errors regions determined by SNIa+CMB+BAO+£/o 
analysis are shown. 



4.2.1 The f(T) PL model 

The power law model [20] assumes the following ansatz of f(T), 

f{T) = a{-T) n . (4.11) 
Here a can be obtained by matching the present matter density f2 m o [20] : 

a = (^"^zt- ( 4 - 12 ) 

Making use of these two equations, Eq.(4.7) can be written as 

E{ z f = n m0 (i + zf + n r0 (i + zf + (i - n m0 - n r0 )E( z ) 2n . (4.13) 

Correspondingly, we also obtain G e // : 

Geff = G r . _ \ rr . (4.14) 

1 l-2n. V H I 

Unlike the ACDM model or the DGP model, this model has two model parameters n 
and f2 m o- In the left panel of Fig.2, we plot the contours of 68.3% and 95.4% confidence levels 
in the {l m Q-n plane, for the f(T)pz model. Constraints from SNLS3, CMB+Hq, BAO+Hq, 
and SNLS3+CMB+BAO+i/o are shown in contours with different colors. At the 68.3% and 
95.4% confidence levels, 

O — fl orko+0.0208 +0.0356 „ _ _ n 9n c:7+0.2490 +0.3663 h _ n 79/ <o+0.0228 +0.0376 ( A i c\ 
"m0 — U.^Ud^_ 0190 _ Q 308) n — U - ZUJ ' -0.3714 -0.7433' " — u - '^°-0.0222 -0.0366' V 1 - 10 ) 

Unlike the DGP model, the different contours given by different observations overlap, showing 
a consistent fit. 

In the right panel of Fig.2, we plot the evolutions of f(z) predicted by the ACDM model 
and the f{T)pL model along with the observational data of growth factor f ^ s . We see that 
the predicted evolutions of f(z) in these two models slightly differ from each other in the la 
CL, especially at the low-redshift region. But the current growth factor data are still not 
powerful enough to distinguish these two models. 
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Figure 3. Left Panel: Probability contours at the 68.3% and 95.4% confidence levels in fl m o-p 
plane, for the f(T)sxp model. Right Panel: The evolutions of f(z) predicted by the ACDM and 
fiT) exp models. The best-fit values and 2a errors regions determined by SNIa+CMB+BAO+f/o 
analysis are shown. 



4.2.2 The f{T) EX p model 

Another popular f(T) model is of exponential form model proposed by Linder [21]. It takes 
the form 



f(T) = cT (l - exp(- PV ¥/T )), (4.16) 

with 

c= - 1 ~° m ° , T = 6Hl (4.17) 

l-(l+p)-p' v ' 

Combining the above two equations with Eqs. (4.7-4. 10), after some tedious calculation, one 
can obtain the E{z) and G e ff of this model, 

E 2 (z) = n m0 (l + zf + n r0 (l + zf + c[l - (1 + pE(z)) eM-pE(z))}, (4.18) 

G ef f = l - T = . (4.19) 

1 _ p(p+l)Pe-V^P(n m0 -l) 



2^/E{z)((p+l)v-l) 

Like the /(T)pl model, this model is also a two-parameter model with parameters p 
and 0, m Q. The contours of the J(T)exp model are shown in the left panel of Fig. 3. As before, 
constraints from different cosmological data are shown in contours with different colors. At 
the 68.3% and 95.4% confidence levels, we obtain 

O — PI 9fi/Lfi+ - 0210 +0-0355 _ n 9 o m +0.4044 +0.8346 x. n 7941 +0.0223 +0.0366 ia r,\ 

s 'm0 — U.^U<±U_o.oi92 -0.0305' P ~ u -^ oui -0.2537 -0.3788' n ~ u - ' z ^ i -0.0216 -0.0354- ^- ZU J 

In the right panel of Fig. 3, we plot the evolutions of f(z) predicted by the ACDM model 
and the f(T) exp model along with the observed data of growth factor f t, s . This figure shows 
that the predicted evolutions f(z) from the expansion history data overlap. Especially, the 
constrained region of f(z) in the ACDM model is a subset of that of the J(T)exp model. 
This means that it will be more difficult to distinguish the cosmic growth history of these 
two models from the growth factor data. 



- 12 - 



4.3 f(R) models 

At last, let us investigate the f(R) model. The basic idea of this model is replacing R by 
R + f(R), yielding the action 

S = TfT~n I d 4 x^ [R + f(R)} + S m + S r , (4.21) 
lb7rG J 

where S m and S r are the actions for the matter content and the radiation content respectively. 
For the background FRW metric, the Ricci scalar can be determined by the Hubble parameter 
and its time derivative, i.e., 

R = 12H 2 + 6H. (4.22) 

Taking variations of the action of f(R) model with respect to the metric in the spatially flat 
FRW universe, one can obtain the modified Friedmann equation: 

H 2 - f R (§ -H 2 ) + f - + H 2 f RR R =^( Pm + Pr). (4.23) 
ob 3 

Here, primes denote derivatives with respect to In a, and f R and f RR are defined by 



The G e ff of the f(R) model is given by [101] 



I fe 2 fRR 



Notice that for the f(R) model, G e ff depends on not only the scale factor but also the 
comoving wavenumber k. As mentioned in [102], the subhorizon approximation cannot be 
satisfied and the non-linear effects are obvious in scale smaller than k = 0.2hMpc~ , while 
the current observations are not so accurate for scale larger than k = 0.01/iilfpc -1 . Therefore, 
for simplicity, we just take k = Q.lhMpc . 

For the f(R) models, we consider two metric forms, proposed by Hu-Sawicki [18] and 
Starobinsky [19], respectively. We will call them J(R)hs an d f(R)st in the following con- 
text. Both of them can satisfy the cosmological and local gravity constraints [104]. In the 
following two subsections, we will discuss in detail the explicit formulas and cosmological 
interpretations of these two models. 

4.3.1 The f(R) H s model 

The model proposed by Hu-Sawicki [18] has the following form of f(R), 

m = R-,R c J^- v (426) 

where \x and n are positive numbers, and R c is the order of the present Ricci scalar Rq. In 
this paper, we take Hu-Sawicki's suggestion of R c = (l m QH^ and further set O m o = 0.25. As 
shown in [103], by using the constraints from the violations of weak and strong equivalence 
principles, Capozziello and Tsujikawa give a bound n > 0.9. In practice [102], n is often 
treated as an integer. For simplicity, we will focus on the case of n = 1. The effects of 
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Figure 4. Left Panel: Probability contours at the 68.3% and 95.4% confidence levels in fi m o-jU 
plane, for the f(R)ns model. Right Panel: the evolutions of f(z) predicted by the ACDM and 
f(R)ns models. The best-fit values and 2a errors regions determined by SNIa+CMB+BAO+£/o 
analysis are shown. In the plot of the f(R) model, k = O.lhMpc^ 1 is chosen to satisfy the subhorizon 
approximation. 

different n will be shown in Appendix B. So actually for the f(R)iis model, we only have 
two free model parameters, i.e., [x and O m o- Correspondingly, the Hubble parameter H(z) 
satisfies the following equation: 




(4.27) 



One can solve this equation numerically to obtain the evolution of H(z). From Eqs.(4.25) 
and (4.26), the G e jf of this f{R) model can also be obtained. 

Now let us discuss the cosmological constraints of the f(R)HS model. In the left panel 
of Fig.4, we plot the contours of 68.3% and 95.4% confidence levels in the fi m o-// plane for 
the /(R)hs model in the case of n = 1. Constraints from SNLS3, CMB+Hq, BAO+Hq, 
and SNLS3+CMB+BAO+-£/o are shown in contours with different colors. At the 68.3% and 
95.4% confidence levels, we obtain 

O _ n 9R/1 o+0.0212 +0.0359 ,, _ -, 7 9Rn o+2.3439 +3.0800 h _ n 71 9n +0.0164 +0.0267 (a no\ 
s ™ — u -^°^°-0.0192 -0.0305' A* — 1 ' • zuuo -3.1935 -3.7323' n ~ u - ' lzu -0.0165 -0.0270' 
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In the right panel of Fig. 4, we plot the evolutions of f(z) predicted by the ACDM model 
and the f(R)ns model along with the observed data of growth factor f b s . These two models 
have similar evolutions of f(z), and it is quite difficult for us to distinguish these two models 
from the current growth factor data. 

4.3.2 The f(R) S t model 

Starobinsky also proposed a famous viable f(R) model [19], in which 



f(R) = R-XR s [l-(l + 



Rl 



(4.29) 



0.25#g. 



where A and n are positive numbers. The same as the /(R)hs model, we choose R s 
As shown in [103], this model also has a bound n > 0.9. So in this work, this model is 
also treated as a two-parameter model with parameter ji and fi m o- Combining the above 
equation with Eqs. (4.23-4.25), the evolution of Hubble parameter H(z) can be obtained by 
numerically solving the following equation: 



1 

+ 6 



R - p,R s 



2\nR 



\2n 



(R/R 



)2n + 1 

-n-1 



+ i 



2XnR, 



+ 



R 2 



R.i 



+ 1 



R 
6 



H 2 ) 



[(2n + l)R 2 




8ttG 



{Pm. + Pr 



(4.30) 



Substituting Eq.(4.29) into Eq.(4.25), one can also obtain the G e // of this f(R) model. 

In the left panel of Fig. 5, we plot the contours of 68.3% and 95.4% confidence levels 
in the Q m o-A plane for the f(R)st model with n = 1. Constraints from SNLS3, CME>+Hq, 
BAO+ffci; an d SNLS3+CMB+BAO-|-i?o are shown in contours with different colors. At the 
68.3% and 95.4% confidence levels, 



(I 



m0 



0.2646 



+0.0215 +0.0363 
-0.0190 -0.0302' 



A = 17.277ll|^2 +3 6433 h = a?122 



+0.0160 +0.0263 
-0.0167 -0.0271 



. (4.31) 



In the right panel of Fig. 5, we plot the evolutions of f{z) predicted by th ACDM model 
and the f(R)st model along with the observational data of growth factor f f, s . Again, we 
obtain similar results: the evolution of the growth factor in the f(R)st model is indistin- 
guishable from that of the ACDM model. 

4.4 Summary 

A brief summary of our results is shown in Table 2, where the models, model parameters, 



number of free model parameters n p , 
ABIC are defined as, 



s, AAIC and ABIC are given. The AAIC and 



AAIC = AlCmodci - AICacdm, 
ABIC = BICmodci — BICacdm- 



(4.32) 
(4.33) 
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Figure 5. Left Panel: Probability contours at the 68.3% and 95.4% confidence levels in fi TO o-A 
plane, for the f(R)st model. Right Panel: The evolutions of f(z) predicted by the ACDM and f(R)st 
models. The best-fit values and 2a errors regions determined by SNIa+CMB+BAO+iJo analysis are 
shown. Again, we choose k = O.lhMpc^ 1 in order to satisfy the subhorizon approximation. 



Table 2. The Summary of Results for five MG models and the ACDM model. 



Models 


Model parameters 


n p 


x 2 ■ 


AAIC 


ABIC 


ACDM 




1 


424.911 








DGP 




1 


470.231 


45.32 


45.32 


f(T)p L 


n, Q m Q 


2 


423.506 


0.595 


4.765 


f(T)EXP 


P, ^mO 


2 


423.275 


0.364 


4.534 


f(R)HS 


^, ^mO 


2 


424.901 


1.990 


5.434 


f(R)st 


A, Q m 


2 


424.903 


1.992 


6.161 



Note: The nuisance parameters h, a, j3 used in the analysis are actually not model parameters with significant 
meanings, so we do not list them in this table. 



The nuisance parameters h, a and (3 are actually not model parameters with significant 
meanings, and hence are not listed in the table. 

To make a comparison, we also list the case of the ACDM model. As shown in the Table 
2, the values of AAIC and the ABIC of DGP model are both quite lager than 6. This means 
that the DGP model is strongly disfavored by the data. Other MG models do not yield to 
remarkable reductions of the Xmin> an d g ive slightly larger AIC and BIC values compared 
with the ACDM model. This indicates that the ACDM model is still more favored by the 
current data. This result is consistent with some previous works [29, 62, 63, 67-69, 105]. 

5 Concluding Remarks 

In this work, we test 5 MG models by using the current cosmological observations. Utiliz- 
ing the observational data of the cosmic expansion history, including the recently released 
SNLS3 type la supernovae sample, the cosmic microwave background anisotropy data from 
the WMAP7 observations, the baryon acoustic oscillation results from the SDSS DR7 and 
the latest Hubble constant measurement utilizing the WFC3 on the HST, we constrain the 
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parameter spaces of these models. Then, by plotting the evolutions of these models' growth 
factor, we further compare the theoretical predictions of these MG models with the current 
growth factor data. It is found that these MG models do not lead to appreciable reductions of 
the Xmim an d give larger AIC and BIC values compared with the ACDM model. In addition, 
based on the current growth factor data, these MG models are difficult to be distinguished 
from the ACDM model, so further growth factor data is needed. 



Acknowledgements 

This work was supported by the NSFC grant No.l0535060/A050207, a NSFC group grant 
No. 10821504. QGH was also supported by the project of Knowledge Innovation Program of 
Chinese Academy of Science and a grant from NSFC (Grant No. 10975167). 



Appendix A. Initial conditions of the matter density perturbation equation 

In this Appendix, we explain the reason of taking initial conditions f(z = 30) = 1 when 
solving the Eq(3.25) numerically. 

Starting from Eq.(3.23) and changing variables from t to a, one can obtain 

d 2 5 | ^,3 | dhiH ^dd A.-KG e ffp m ^_^ 
da 2 a da da a 2 H 2 

To solve this equation numerically, we take the initial condition at high-z era, e.g. z = 30 
here, since G e jf = G is satisfied precisely at high-2 era for all models considered in this 
paper. 

In addition, in the high-z regime, the universe is at the matter-dominated stage, thus 
we have the Friedmann equation: 



8vrG 



Then, Eq.(5.1) becomes 



H z = —Pm- (5.2) 



^ S 3 d6 3 r n te oN 

^ + 2^-2^ = °- (5 ' 3) 



Assuming the solution of the above equation takes the form 5 = 5 cons ta m , and then 
substituting this form into Eq.(5.3), one finally gets the function's general solution 

5 = C ia + C 2 a- 3/2 . (5.4) 

Since 5 is quite small in high-z regime, the acceptable solution is 

5 = da. (5.5) 

Thus, we finally obtain the initial condition of Eq.(3.25) 

din 5 , 

f = -77— = 1 5.6 
a ma 

at redshift z = 30. 
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Figure 6. The evolutions of f(z) predicted by the ACDM and /(R)hs models in the cases 
of n = 1 and n = 2 respectively. The best-fit values and 2er errors regions determined by 
SNIa+CMB+BAO+l/o analysis are shown. Note that, k = O.lhMpc' 1 was chosen in order to 
satisfy the subhorizon approximation. 

Appendix B. The effects of n in the two f(R) models 

In this Appendix, we show that the effects of n on the cosmological interpretations of the 
f(R) models are rather small. 

Let us consider the f(R)ns model. The best-fit values of the cosmological parameters 
O m o and h, their 68.3% and 95.4% confidence region, and the corresponding x m in can be 
obtained from the joint analysis of the SNLS3+CMB+BAO+-ffo data. When we take n = 1, 
we find the following results, 

O — 9R/L8+ - 0212 +0-0359 h _ n 71 9n +0.0164 +0.0267 2 _ A0A Qm /c 7 \ 

1! ">0 - U.2b48_ _ 0192 _ 0i0305 , h - U.712U_ 0165 _ ao27 o, Xmin - 424.yUl, (b.l) 

while n = 2 yields to 

O — fl 9fiA8+ a ° 214 +0-0357 h _ n 7191 +0.0160 +0.0265 ^2 _ Ar>A Q1 o /c o\ 

iZ m0 - U.2b48_ 0190 _ . 0305 , h - U.a21_ 00165 „ .0268' Xmin ~ 424.yid. (t>.8J 

Clearly, the above two sets of results are very close to each other. 

In Fig. 6, we also plot the evolutions of f(z) predicted by the ACDM model and the 
f(R)iis model with n = 1 and n = 2, along with the observed data of growth factor / {, s . 
Also, it can be seen that the results of n = 1 and n = 2 are similar to each other. 

In all, it is clear that different n give quite similar results of the constraints on f2 m o 
and h, the XminS, and the evolutions of f(z). That is, the effects of n on the cosmological 
interpretations of the /(R)hs model are rather small. In the f(R)st model, the result is 
similar. Therefore, as mentioned above, it is unnecessary for us to treat n as a free model 
parameter for the two f(R) models. 

References 

[f] A. G. Ricss et al., AJ. 116, 1009 (1998); S. Perlmutter et al., ApJ. 517, 565 (1999). 



- 18 - 



[2] P. J. E. Peebles and B. Ratra, Rev. Mod. Phys. 75, 559 (2003); T. Padmanabhan, Phys. Rept. 
380, 235 (2003); E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. Mod. Phys. D 15, 1753 
(2006); A. Albrecht, et al., astro-ph/0609591; J. Frieman, M. Turner and D. Huterer, Ann. Rev. 
Astron. Astrophys 46, 385 (2008); S. Tsujikawa, arXiv:1004.1493; M. Li et al, arXiv:1103.5870. 

[3] B. Ratra and P.J.E. Peebles, Phys. Rev. D 37, 3406 (1988); P. J. E. Peebles and B.Ratra, ApJ 
325, L17 (1988); C. Wetterich, Nucl. Phys. B302, 668 (1988); I. Zlatev, L. Wang and P. J. 
Steinhardt, Phys. Rev. Lett. 82, 896 (1999). 

[4] R. R. Caldwell, Phys. Lett. B 545, 23 (2002); S. M. Carroll, M. Hoffman and M. Trodden, Phys. 
Rev. D68, 023509 (2003). 

[5] C. Armendariz-Picon, T. Damour and V. Mukhanov, Phys. Lett. B 458, 209 (1999); C. 
Armendariz-Picon, V. Mukhanov and P. J. Steinhardt, Phys. Rev. D 63, 103510 (2001); T. 
Chiba, T. Okabe and M. Yamaguchi, Phys. Rev. D62, 023511 (2000). 

[6] A. Y. Kamcnshchik, U. Moschella and V. Pasquier, Phys. Lett. B 511, 265 (2001); M. C. Bento, 
O. Bcrtolami and A .A. Sen, Phys. Rev. D 66, 043507 (2002); X. Zhang, F. Q. Wu and J. Zhang, 
JCAP 0601 003 (2006); S. Li, Y. G. Ma and Y. Chen, Int. J. Mod. Phys. D 18 1785 (2009). 

[7] T. Padmanabhan, Phys. Rev. D66, 021301 (2002); J. S. Bagla, H. K. Jassal, and T. 
Padmanabhan, Phys. Rev. D 67, 063504 (2003). 

[8] M. Li, Phys. Lett. B 603, 1 (2004); Q. G. Huang and M. Li, JCAP 0408, 013 (2004); Q. G. 
Huang and Y. G. Gong, JCAP 0408, 006 (2004); Q. G. Huang and M. Li, JCAP 0503, 001 

(2005) ; X. Zhang and F. Q. Wu, Phys. Rev. D 72, 043524 (2005); B. Wang, E. Abdalla and R. K. 
Su, Phys. Lett. B 611, 21 (2005); B. Wang, C. Y. Lin and E. Abdalla, Phys. Lett. B 637, 357 

(2006) ; J. Zhang, X. Zhang and H. Y. Liu, Eur. Phys. J. C 52, 693 (2007); C. J. Feng, Phys. Lett. 
B 633, 367 (2008); Y. Z. Ma, Y. Gong and X. L. Chen, Eur. Phys. J. C 60, 303 (2009); M. Li et 
al, Commun. Thcor. Phys. 51, 181 (2009); M. Li et al, JCAP 0906, 036 (2009); M. Li et al, 
JCAP 0912, 014 (2009); X. Zhang, Phys. Lett. B 683, 81 (2010). 

[9] H. Wei and R. G. Cai, Phys. Lett. B 655, 1 (2007); R. G. Cai, Phys. Lett. B 657, 228 (2007); H. 
Wei and R. G. Cai, Phys. Lett. B 660 113 (2008); H. Wei and R. G. Cai, Phys. Lett. B 663, 1 
(2008); J. Zhang, X. Zhang and H. Liu, Eur. Phys. J. C 54, 303 (2008); J. P. Wu, D.Z. Ma and Y. 
Ling, Phys. Lett. B 663, 152 (2008). 

[10] S. Nojiri and S. D. Odintsov, Gen. Rel. Grav. 38, 1285 (2006); C. Gao, F. Wu, X. Chen and Y. 
G. Shcn, Phys. Rev. D 79, 043511 (2009); C. J. Feng, Phys. Lett. B 670, 231 (2008); C. J. Feng, 
Phys. Lett. B 672, 94 (2009); L. N. Granda and A. Oliveros, Phys. Lett. B 669, 275 (2008); X. 
Zhang, Phys. Rev. D 79, 103509 (2009); C. J. Feng and X. Zhang, Phys. Lett. B 680, 399 (2009). 

[11] H. Wei, R. G. Cai, and D. F. Zcng, Class. Quant. Grav. 22, 3189 (2005); H. Wei, and R. G. 
Cai, Phys. Rev. D72, 123507 (2005). 

[12] Y. Zhang, T. Y. Xia, and W. Zhao, Class. Quant. Grav. 24, 3309 (2007); T. Y. Xia and Y. 
Zhang, Phys. Lett. B 656, 19 (2007); S. Wang, Y. Zhang and T. Y. Xia, JCAP 10, 037 (2008); S. 
Wang and Y. Zhang, Phys. Lett. B 669, 201 (2008). 

[13] V. K. Oncmli and R. P. Woodard, Class. Quant. Grav. 19, 4607 (2002); V. K. Oncmli and R. 
P. Woodard, Phys. Rev. D 70, 107301 (2004); E. O. Kahya and V. K. Oncmli, Phys. Rev. D 76, 
043512 (2007). 

[14] S. Weinberg, Rev. Mod. Phys. 61, 1 (1989); S. M. Carroll, W. H. Press and E. L. Turner, Ann. 
Rev. Astron. Astrophys. 30, 499 (1992); V. Sahni and A. Starobinsky, Int. J. Mod. Phys. D 9, 373 
(2000); S. M. Carroll, Living Rev. Rel. 4, 1 (2001). 

[15] S. Nojiri and S. D. Odintsov, Phys. Rept. 505, 59 (2011). 

[16] T. Clifton et al, arXiv: 1106. 2476. 

[17] G. R. Dvali, G. Gabadadze, and M. Porrati, Phys. Lett. B 485, 208 (2000). 



- 19 - 



[18] W. Hu and I. Sawicki, Phys. Rev. D 76, 064004 (2007). 
[19] A. A. Starobinsky, J. Exp. Thcor. Phys. Lett. 86, 157 (2007). 
[20] G. R. Bengochea and R. Ferraro, Phys. Rev. D 79, 124019 (2009). 
[21] E. V. Lindcr, Phys. Rev. D 81, 127301 (2010). 

[22] M. Milgrom, ApJ. 270, 365 (1983); J.Bckenstein and M. Milgrom, ApJ. 286, 7 (1984); J. D. 
Bekenstein, in Proceedings of the Sixth Marcel Grossman Meeting on General Relativity, H. Sato 
and T. Nakamura, eds (World Scientific, Singapore 1992). 

[23] J. D. Bekenstein, Phys. Rev. D 70, 083509 (2004). 

[24] C. Brans and R. H. Dickc, Phys. Rev. 124, 925 (1961); L. Amcndola, Phys. Rev. D 60, 043501 
(1999). 

[25] B. Zwicbach, Pys. Lett. B 156, 315 (1985). 

[26] P. Horava, Phys. Rev. D 79, 084008 (2009); E.N. Saridakis, Eur. Phys. J. C. 67, 229 (2010). 

[27] S. Nojiri and S.D. Odintsov, Int. J. Geom. Mcth. Mod. Phys. 4 115 (2007); T.P. Sotiriou and 
V. Faraoni, Rev. Mod. Phys. 82, 451 (2010); A.D. Felice and S. Tsujikawa, Living. Rev. Rel. 13, 
3 (2010); 

[28] R.K. Sachs and A.M. Wolfe, ApJ. 147, 73 (1967). 
[29] D. Rubin et ai, ApJ. 695, 391 (2009). 

[30] Zong-Kuan Guo, Zong-Hong Zhu, J.S. Alcaniz and Yuan-Zhong Zhang, Astrophys. J. 646, 1 
(2006). 

[31] Matteo Martinelli, Alessandro Melchiorri and Luca Amcndola, Phys. Rev. D 79, 123516 (2009). 
[32] Puxun Wu and Hongwci Yu, Phys. Lett. B 693, 415-420 (2010). 

[33] Amna Ali, Radouane Gannouji, M. Sami and Anjan A. Sen, Phys. Rev. D 81, 104029(2010). 

[34] F.C. Carvalho, E.M. Santos, J.S. Alcaniz and J. Santos, JCAP 0809, 008(2008). 

[35] Abha Dcv, D. Jain, S. Jhingan, S. Nojiri, M. Sami and I. Thongkool, Phys. Rev. D 78, 083515 
(2008). 

[36] J. Santos, J.S. Alcaniz, F.C. Carvalho and N. Pires, Phys .Lett. B 669, 14 (2008). 

[37] Louis Yang, Chung-Chi Lee, Ling- Wei Luo and Chao-Qiang Geng, Phys. Rev. D 82, 103515 
(2010). 

[38] Koichi Hirano and Zen Komiya, Int. J. Mod. Phys. D 20, 1 (2011). 

[39] Hao-Yi Wan, Zc-Long Yi and Tong-Jic Zhang, Phys. Lett. B 651, 352 (2007). 

[40] A. Luc, R. Scoccimarro and G. Starkman, Phys. Rev. D 69, 124015 (2004). 

[41] V. Acquaviva, A. Hajian, D.N. Spergel and S. Das, Phys. Rev. D 78, 043514 (2008). 

[42] Y.G. Gong, Phys. Rev. D 78, 123010 (2008). 

[43] D. Polarski and R. Gannouji, Phys. Lett. B 660, 439 (2008). 

[44] E.V. Lindcr, Phys. Rev. D 72, 043529 (2005). 

[45] K. Koyama and R. Maartens, JCAP 0601, 016 (2006). 

[46] T. Koivisto and D.F. Mota, Phys. Rev. D 73, 083502 (2006). 

[47] S. Daniel, R. Caldwell, A. Cooray and A. Melchiorri, Phys. Rev. D 77, 103513 (2008). 

[48] L. Knox, Y.-S. Song and J.A. Tyson, Phys. Rev. D 74, 023512 (2006). 

[49] M. Ishak, A. Upadhye and D.N. Spergel, Phys. Rev. D 74, 043513 (2006). 



- 20 - 



[50] I. Laszlo and R. Bean, Phys. Rev. D 77, 024048 (2008). 

[51] B. Jain and P. Zhang, Phys. Rev. D 78, 063503 (2008). 

[52] W. Hu and I. Sawicki, Phys. Rev. D 76, 104043 (2007). 

[53] A. Conley et ai, ApJS. 192, 1 (2011). 

[54] J. Guy et ai, arXiv:1010.4743. 

[55] M. Hamuy et ai, A3. 112, 2408 (1996); A. G. Ricss et at, AJ. 117, 707 (1999); S. Jha et ai, 
AJ. 131, 527 (2006); C. Contreras et ai, AJ. 139, 519 (2010). 

[56] M. Hicken, et ai, ApJ. 700, 1097 (2009); M. Hicken, et ai, ApJ. 700, 331 (2009). 

[57] J. A. Holtzman et at, A3. 136, 2306 (2008). 

[58] A. G. Riess et at, ApJ. 659, 98 (2007). 

[59] E. Komatsu et at, ApJS. 192, 18 (2011). 

[60] W. J. Percival et ai, MNRAS 401, 2148 (2010). 

[61] A. G. Riess et at, ApJ. 730, 119 (2011). 

[62] H. Akaikc, IEEE Trans. Automatic Control 19, 716 (1974). 

[63] G. Schwarz, Ann. Stat., 6, 461 (1978). 

[64] A. A. Starobinsky, JETP Lett. 68, 757 (1998). 

[65] A. Lewis and S. Bridle, Phys. Rev. D 66, 103511 (2002). 

[66] A. R. Liddlc, Mon. Not. Roy. Astron. Soc. 351, L49 (2004). 

[67] W. Godlowski and M. Szydlowski, Phys. Rev. Lett. B 623, 10 (2005). 

[68] M. Biesiada, JCAP 0702, 003 (2007). 

[69] J. Magueijo and R. D. Sorkin, MNRAS 377, L39 (2007). 

[70] M. Szydlowski and W.Godlowski, Phys. Lett. B 633, 427 (2006); M. Szydlowski, A. Kurek and 
A. Krawiec, Phys. Lett. B 642, 171 (2006); P. Mukherjee et ai, Mon. Not. Roy. Astron. Soc. 369, 
1725 (2006); 

[71] A. R. Liddle, Mon. Not. Roy. Astron. Soc. Lett. 377, L74 (2007). 

[72] M. Sullivan et at, ApJ. 737, 102 (2011). 

[73] X.D. Li et at, JCAP 1107, 011 (2011). 

[74] Y. Wang, C.H. Chuang and P. Mukherjee, arXiv:1109.3172. 

[75] Y.G. Gong, Q. Gao and Z. H. Zhu, arXiv:1110.6535. 

[76] ZX. Li, P.X. Wu and H.W. Yu, arXiv:1109.6125. 

[77] https://tspace.library.utoronto.ca/handle/1807/24512 

[78] M. Sullivan et at, MNRAS 406, 782 (2010). 

[79] W. Hu and N. Sugiyama, ApJ 471, 542 (1996). 

[80] J. R. Bond, G. Efstathiou and M. Tegmark, Mon. Not. R. Astron. Soc 291, L33 (1997). 

[81] D. J. Eisenstein et ai, ApJ 633, 560 (2005). 

[82] D. J. Eisenstein and W. Hu, ApJ. 496, 605 (1998). 

[83] W. L. Frccdman and B. F. Madore, arXiv:1004.1856. 

[84] W. Hu, ASP Conf. Scr. 339, 215 (2005). 

[85] B. Boisseau, G. Esposito-Farse, D. Polarski, A. A. Starobinsky, Phys. Rev. Lett. 85, 2236 (2000). 



- 21 - 



[86] C. Di Porto and L. Amendola, Phys. Rev. D 77, 083508 (2008). 

[87] S. Nesscris and L. Pcrivolaropoulos, Phys. Rev. D 77, 023504 (2008). 

[88] L. Guzzo et ai, Nature 451, 541 (2008). 

[89] M. Colless et ai, Mont. Not. R. Astron. Soc. 328, 1039 (2001). 

[90] M. Tcgmark el ai, Phys. Rev. D 74, 123507 (2006). 

[91] N.P. Ross et ai, Mont. Not. R. Astron. Soc. 381, 573 (2007). 

[92] J. da Angela et ai, Mont. Not. R. Astron. Soc. 383, 565 (2008). 

[93] P. McDonald et at, Astrophys. J. 635, 761 (2005). 

[94] M. Vicl, M.G. Haehnelt and V. Springel, Mont. Not. R. Astron. Soc. 354, 684 (2004). 

[95] M. Vicl, M.G. Haehnelt and V. Springel, Mont. Not. R. Astron. Soc. 365, 231 (2006). 

[96] C. Deffayet, Phys. Lett. B 502, 199 (2001); C. Deffayet, G. R. Dvali and G. Gabadadze, Phys. 
Rev. D 65, 044023 (2002). 

[97] A. Lue, R. Scoccimarro and G. Starkman, Phys. Rev. D 69, 044005 (2004). 

[98] W. J. Fang, et ai, Phys. Rev. D 78, 103509 (2008); M. Fairbairn and A. Goodbar, Phys. Lett. 
B 642, 432 (2006); R. Maartens and E. Majerotto, Phys. Rev. D 74, 023004 (2006); U. Alam and 
V. Sahni, Phys. Rev. D 73, 084024 (2006); Y. S. Song, I. Sawicki, and W. Hu, Phys. Rev. D 75, 
064003 (2007); F. Schmidt, Phys. Rev. D 80, 043001 (2009). 

[99] G.R. Bengochea and R. Ferraro, Phys. Rev. D 79, 124019 (2009). 

[100] Rui Zheng, Qing-Guo Huang, JCAP 1103, 002 (2011). 

[101] Shinji Tsujikawa, Phys. Rev. D 76, 023514 (2007). 

[102] X. Fu, P. Wu, H. Yu, Eur. Phys. J. C 68, 271 (2010). 

[103] S. Capozzicllo and S. Tsujikawa, Phys. Rev. D 77, 107501 (2008). 

[104] Shinji Tsujikawa, Phys. Rev. D 77, 023507 (2008). 

[105] T. M. Davis et ai, ApJ. 666, 716 (2007); J. Sollcrman et at, ApJ. 703, 1374 (2009); M. Li, X. 
D. Li and X. Zhang, Sci. China Phys. Mech. Astron. 53, 1631 (2010); H. Wei, JCAP 1008, 020 
(2010). 



-22- 



